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ABSTRACT 

The direct problem in transient heat conduction requires the 
determination of conditions at an interior location when conditions 
are known at the boundaries of a solid. Conversely, the inverse 
problem requires the determination of conditions at the boundaries 
of a solid when conditions are known at an interior location.  Con- 
sequently special methods are required in the solution of the inverse 
problem. A new method, numerical inversion of the Laplace transform, 
is used to solve this complex problem. Application of this numerical 

код die abs 

technique to the seni-infinitdseemepeetm с" cylinder, and sphere is 
made, and the accuracy of solution is discussed. This method of 
solution provides the engineer with a simple, powerful tool that can 


be used in the determination of heat transfer phenomena in a solid. 


Thesis by: Terrill Jay Wendt, entitled "Application of Numerical Inversion 


of the LaPlace Transform to the Inverse Problem in Transient Heat Conduction". 
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Ts INTRODUCTION 

Traditionally, heat-transfer problems are concerned with the 
determination of conditions at an interior location when conditions 
are known at the boundaries of a solid. Such problems are known as 
direct problems, and their solutions are obtained through straight- 
forward classical techniques M In many applícations (i.e., 
re-entry thermodynamics, rocket motors, quenching), it is desired 
to determine the conditions at the boundaries of a solid when 
conditions are known at an interior location. Problems of this type 
are known as inverse problems and require special methods of solution. 

Many solutions to the inverse problem of heat conduction exist 
in the literature. Stolz [19] uses numerical inversion of a specific 
direct problem to solve for the heat flux at the surface of a solid. 
He notes that there is appreciable damping or lag in time of the heat 
flux from the interior location to the surface. This damping resulted 
in oscillations of the approximate solution with resultant errors. 
Burggraf [4] obtains an "exact" solution to the inverse problem in 
heat conduction with internal heat generation. However Beck [4] in 
commenting on this solution notes that this "exact" method must 
utilize inexact measured temperature data for which accurate time 
derivatives (necessary for this "exact" solution) are very difficult 
to obtain. In commenting on the inaccuracies involved in using 
measured temperature data, Frank [12] suggests using a least squares 
fit to the data available in order to obtain an essentially smooth 
temperature function for use in the solution of the inverse problem. 


Sparrow, et. al. [18] use a modification of Laplace transform inversion 


* 
Numbers in brackets refer to references in the bibliography. 
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by applying a convergence factor to the temperature function in the 
S plane. Deverall and Channapragada [11] utilize a convergence factor 
to obtain a solution for the heat flux in an inverse problem. 

It is intended here to present a simple straightforward method 
of obtaining solutions to inverse problems of practical interest. The 
numerical technique involved in inverting the Laplace transform (used 
in solving the inverse problems) was developed by Bellman and Kalaba 
[2 and 3]. This method, based on orthogonal function theory and 
Gaussian quadrature, is developed and tested for accuracy (Tables I 
and II). This numerical technique is applied to the semi-infinite 
solid, the "long" cylinder and the sphere. The accuracy of solution 
is compared graphically with the known exact solutions. The FORTRAN IV 
programs used in the solutions are listed in Appendix A for use by 


the reader. 


ша 


TT. NUMERICAL TECHNIQUE 
INTRODUCTION 

The numerical technique involved was developed by Bellman and 
Kalaba [2] and involves the application of orthogonal function theory 
to the inversion of the Laplace Transform as opposed to using the 
Convolution Theorem which involves integration in the complex plane. 

It is important to note that the method developed in this 
section can be applied to other problems in physics, biology, eco- 


nomics, chemotherapy and radiation [2]. 


THE LAPLACE TRANSFORM 
We will first consider some of the properties of the Laplace 
transform and limitations involved in applying the inverse. 
Let f(t) be any sectionally continuous function defined for 
t > 0. Then the Laplace transform of f(t) is defined as 


со 


F(s) = L{f(t)} = | cali f(t)dt (1) 


O 


where the real part of s is large enough to make the integral (1) 
convergent. 

The Laplacian Operator L is a linear operator and possesses 
many desirable properties, mathematically and computationally. Some 
of these properties are illustrated in the following theorems which 
are included here for completeness: 


THEOREM I: ШЕ L(f (t)) - F (s) and L(f (t)} = F (s) 


then L{f (t) + f (t)} = F (s) + F (s). 
1 2 1 2 


THEOREM II: If L{f(t)} = F(s) 


then 1850) = sF(s) - £(0). 
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THEOREM III: 


THEOREM IV: 


THEOREM V: 


THEOREM VI: 


THEOREM VII: 


THEOREM VIII: 


ТЕ ЕСЕ) а (в) 


N N 

кесе Цецо) _ I 
3x dx 
Ж x 


Tf L{f(t)} = F(s) 


arc мам! 
then {| көзе) a F(s). 
о 


If L{f(t)} = F(s) 


then L{f(at)} = rÊ). 


Initial Value Theorem. 

ТЕ та Ее) е Е(в) 
then 

E әк (ЕЕ Linis 


Final Value Theorem. 


If LE OE ECs) 


then 
Limit Limit 
о sF(s) = mus f(t). 


Superposition Theorem. 
TE bitte) = (57 and 


then 


t 
4 Osca] = 4 
О 


T EUR 


Lig(t)} = G(s) 


t 
osea 


О 


"^ F(s)G(s) 
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INSTABILITY OF THE INVERSE OF THE LAPLACE TRANSFORM 
The instability of the inverse of the Laplace transform can 
be illustrated by considering the function 
f(t) = sin Ot. 
The Laplace transform of f(t) is then 


L{sin wt} = (2) 


EU 
52-2 
Now it is known that sin ot is bounded in amplitude by (1,-1); but as 
Q is increased, sin wt oscillates correspondingly more rapidly. 
However the Laplace transform of sin wt approaches zero as w is 


increased. 


Limit 0) = 0 


52-2 
(0-00 


In other words small changes in s can cause arbitrarily large 
changes in t with respect to the period of oscillation. Consequently 
we cannot sift or filter out rapid oscillations by using values of 
F(s). Such difficulties as mentioned above can be avoided only if 
f(t) is an essentially smooth function, which is an implicit assump- 


tion in the following development. 


NUMERICAL QUADRATURE 
There are several methods that can be used to approximate a 
definite integral in terms of a finite sum. The method employed here 
is based on Gaussian Quadrature and is used because of convergence 


properties and exactness of solution [15]. Consider the integral 


approximation 
b M 
| w(x) £(x)dx = за А Е(х.). (3) 
1 1 
а 1=1 
By minimizing the function I 
M 
wx)£(x) - > Е (4) 


15 


we can assume that the summation involved in equation (3) is exact 
for a prescribed accuracy, €, for any polynomial of degree less than 
or equal to M. The х, 5 and the А !8 are to be determined so that 
equation (3) is exact for a polynomial f(x) of high degree. Here 
the х,'8 are called the roots (or nodes), the A,'s are called the 
coefficients (or weights), and w(x) is called the weighting function 
of the function f(x). 

With this background let us examine the numerical approxima- 


tion to an integral which will be used later to invert the Laplace 


transform. Let 


M 
T(x) = T] (x-x1) , (5) 
i=1 
and let Р(х) be a polynomial of degree N. Then we have 
РО) 7 TGOQy, t Py- (6) 
Since Р(х.) = Р, 104) TS l 205 М 


b 


b b 
then | и(к)Р „Содх = | w(x)P_ 4 (x) dx + | wGom(oQ, , (x)dx. (7) 
а а 


а 


Now let us choose the x, 'S such that 


b 
| ыб) п(к)0, у (х)дх = 0, (8) 
а 


for all polynomials 0, м09. Then from equation (3) we have 


b 
| Угус = E АР. (к.). (9) 


i=1 


If equation (8) is to be valid for all polynomials Qs. M for 


N > M, it must be valid гог 
D K 
| w(x)T(x)x dx = 0, к = 1,2,...,N-M (10) 
а 
Let us now expand T(x) into a polynomial 
i 


T(x) = ES Сх e. (11) 
i=o 
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so that equation (10) becomes 


b M | 
11 К š 
K a с) dx = 0, і = 0,1,2,...,N-M (12) 


which is a system of M equations in the с. 8 if N = 2M-1. 

The system of equations in (12) has a unique solution if and 
only if the matrix of coefficients is nonsingular. 

The determinant of this matrix is called a Vandermonde deter- 


minant and can be easily evaluated, since w(x) > 0 in the inverval 


(a,b), as 
1 x x ы 
O O 
N 
1 X, š X, 
= (х.-х.). (13) 
' а pj ^ J 


1 N 
XN е е e XN | 
Since the х 8 are real and distinct points (roots), the determinant 


in (13) does not vanish and the b,'s can be uniquely determined. 


i 


It can be shown that the zeros of T(x), determined by its 
coefficients (c1), are real and in the interval (a,b). It can also 
be shown that Qu X. C9 is a polynomial of degrees M-1, where N = 2M-1. 

Now let T, (x) and n designate the functions Tx) deter- 


mined by their coefficients for M = i and M = j respectively. Then 


b 
: | сл сс = 0 i = j (14) 
a 


If i <.3, then T, (x) is a polynomial of degree less than j; and the 
result follows from equation (8), which has been verified. From 
equation (14), T, (x) and T are said to be orthogonal on the 


interval (a,b) with respect to the weight function w(x). 
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For each value of M, the х 8 in equation (3) are taken as 
the zeros of TOO as determined by its coefficients. The A,'s in 
equation (3) can then be determined by the method of undetermined 
coefficients. It can be shown that all of the А в are positive. 

Therefore for a given weight function w(x) and interval of 
integration (a,b), the zeros of the polynomials Ty 08) » as determined 
by their coefficients (с), need be computed only once. Likewise 


° ! 
the corresponding А, s need be computed only once. 


GAUSS-LEGENDRE QUADRATURE 

Since we have been dealing with any orthogonal function mu) , 
let us indicate the specific orthogonal functions used in the fol- 
lowing development. The Legendre polynomials are the unique Gauss 
quadrature polynomials [15] in the interval (-1,1) (or any finite 
interval) with weighting function 

VOCE (15) 

Let us denote the Legendre polynomials, in accordance with 


common practice, by the notation PAG where 


2N)! 
рб) = та) NAN > (16) 


and the interval of expansion is (-1,1). These polynomials are defined 
by the Rodrigue's formula 


D E a^ ei (17) 
N 2NN! arN i 


From equation (14) it follows that 





EN Св — 5 MTN. (18) 
-1 
The recursion relations for the Legendre polynomials are then 
es (241); EN = Gr en oe CZ 
r^-1 ' 
and Pau 1 EE P(r) + Del LN (r). (20) 
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Consequently the least squares polynomial approximation 
(Gaussian Quadrature) to f(x) over (-1,1) with respect to constant 


weighting function is 





N 
f(r) = > b P (r), (21) 
K-0 
1 
2K+1 
where Ь a EX (22) 
Therefore the coefficients (or weights) А, аге 
А | ш d (23 
aa) =e See. > 
i т г.) P (r) 


where r.'s are the roots of the Legendre polynomials Ры (к). 


SHIFTED LEGENDRE POLYNOMIALS 

We have now examined a specific orthogonal polynomial, the 
Legendre polynomial. However, it is most convenient in the following 
development to consider the shifted Legendre polynomials. The reason 
for employing these polynomials is twofold. First, the interval of 
expansion (0,1) of the shifted Legendre polynomials corresponds to 
that involved in obtaining an explicit inverse matrix. Second, the 
coefficients of the shifted Legendre polynomials are all integers, 
which is an important property computationally. 

Since the roots and weights have been tabulated for the 
interval (-1,1) [19], it is a simple matter to obtain the roots and 
weights of the shifted Legendre polynomials in the interval (0,1). 


Let х 
PNG = 2 (24) 


* 
where PN (x) are the shifted Legendre polynomials. Consequently the 


roots of the shifted Legendre polynomials are 





Iera 
X, = › (25) 
i 2 
and the weights are 
х Ai 
А, = 75 · (26) 


T3 


NUMERICAL INVERSION OF THE LAPLACE TRANSFORM 
We have now laid a proper foundation for obtaining an explicit 
inverse matrix to be used in inverting the Laplace transform. Let us 
consider the approximate inversion of the Laplace transform through 
numerical techniques. 


Let M 
-st 
F(s) = | e u(t)dt. (27) 
O 
Now we make a change of the variable of integration to 
х= е" (28) 
therefore t-- log x (29) 


Substituting equation (29) into equation (27) we obtain 


1 
| zonis x)dx = F(s). (30) 
О 
Now let a(x) = u(-log x) (31) 


and then equation (30) becomes 


1 — 
| NR dees Cae (32) 
О 
If we make g(x) = а(хух,871, then using the quadrature formula 
1 N a 
| g(x)dx = > A, g(x,), (33) 
с 1 T 
i=l 
О 
we obtain 
S Ж 5-1 
>_ А, X. а(х.) = F(s). (34) 


i 
Letting s = 1,2,...,N, we obtain a linear system of N equations in 


the N unknowns, а(х.), ism 12 NS 
N 


* 
T9 A. x a(x,) = Е(«+1) K 2 0,1,...,N-1 (35) 
2 1 1 
1=1 
It is our intention here to obtain 
N-1 
а(х.) = > аР), dom S N. (36) 


where the E are constants independent of F(K-*1), determined once 


and for all for each N, and tabulated. 
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Let us now make another change of variables such that 


* 
У = A, ox), (37) 
and B s Es (38) 
Equation (35) then becomes 
N 
K pe 
> хуу; = а, к 0, ] ПЕШ (39) 
i=l 


Multiplying the po equation by a parameter 4, and adding the 


corresponding terms in all the equations, we get 


N -1 N-1 
N K 
> y, > 4,х, = > ad (40) 
i=l K-o С-о 
Now setting 
N-1 = 
f (x) = > Hes 5 (51) 
K=O 
N N-1 
we have 2 y, fog = >а, К (42) 


where f(x) is a polynomial of degree №1, to be chosen in some 
convenient manner. 
We want to determine the Y, 9» and thus we require that 


LONE Е С) such that 


o К 

an orthogonality condition. If 
f (x) -> dea X^ o (44) 

equation (42) reduces to 


YS vc (45) 


He þe 
lI + 
Ca. Ca. 


° (43) 


ca 
^ 
il 
о o 


where eq? K = 0,1,...,N-l, are determined by the orthogonality 
condition in equation (43). 
Since the interval of integration involved in the preceding 


development is (0,1), we will use the shifted Legendre polynomials 


in approximating the integral in equation (32). Because the x, 'S are 
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х 
the zeros of the shifted Legendre polynomial Ры (x), we can use the 


Lagrange interpolation formula 


И) 
i (x) = te : (46) 


The desired 4.4 5 in equation (45) are the coefficients in 
this polynomial of degree N-1. Repeating this procedure for each j 


yields the desired inverse matrix (q .). 
KJ 


INTERLACING OF THE ROOTS 
It should be noted that the roots (х.) for different N inter- 
lace (Appendix B). Consequently the values of time 
t, = -log х; (47) 
intertwine, and the results for one N may be combined with those of 
another N to obtain both a check on the results and an enlarged set 


of t values for extrapolation purposes. 


CHANGE OF TIME SCALE 
The method of inversion described for the Laplace transform 
provides approximate values for the function u(t) at discrete values 


of t. For general N, we obtain u(t.) at the points t. = -log X. where 


* 
(x). 


x, are the roots (or zeros) of the shifted Legendre polynomials, P 


i N 


It is obvious that the inverse transform is thus determined for 
a very small range of t values. Since increasing N does not signifi- 
cantly increase the range of t, we must resort to another method of 
time scaling. Recall that if 


со 
F(s) = L{u(t)} = | e ?*u(t)dt, (48) 
О 
then from Theorem V, it follows that 


L{u(at) } - =F). (49) 
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Therefore using the approximate relation 


ч 5-1 1 S 
> Ах. u(-a log x.) Зла (=) i (50) 
і-1 

we can determine u(t) for values of A - -a log X; 


ALTERNATE METHODS 
Although the method of inverting the Laplace transform pre- 
sented involves the use of the shifted Legendre polynomials, there 
are several other methods that can be used in obtaining an explicit 
inverse matrix. The first is the Method of Papoulis which utilizes 
the Legendre polynomials and the second is based on a different set 
of orthogonal polynomials, the Laguerre polynomials. 
Papoulis Method 
The Papoulis method of inversion involves the expansion of f(t) 
into a series of Legendre polynomials and uses the classical interval 
of orthogonality (-1,1). It is clear from this method how a particular 
region of the time scale may be examined by changing the scale factor 
ү [8]. 


Let us make the transformation 


x = ae Y > 0 (51) 
where f(t) becomes 
r ton s Е ае s (52) 
Making these substitutions into equation (1) we get 
F(s) - hte f (x) dx : (53) 
O 


By expressing s in terms of y as s = (2kK+1)y, equation (53) becomes 


1 
ҮЕ(2кү+ү) = | NS 


О 


"f(x)dx. (54) 
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Defining f (x) as an even function (i.e., f (x) - f(-x)), the function 


f (x) can be expanded into a series of even Legendre polynomials as 


f(x) = > сұр, (е). (55) 
N20 
EXEC. : : -Yt 

Using the fact that Р „(е ) is an even polynomial in e 
of degree 2N, Papoulis shows that the transform of P, (е ТО) Ш5 МО 
the form 

- _ (s-y) (s-3y) ...[s-(2N-1) y] 
F(s) Pon's) S(st2y) ... (st2Ny) Ca 


Replacing s first by y then 3y, 5y,...(2N+1)y ... in equation (56), 


gives the following system of equations to evaluate Сү: 


YEGY) СС 
Т 2 
УЕ (Зу) = 3 2- + 15 Ci 
(505 
| Co 2NCi 2N(2N-2) ...2CN 
FONYTY) 7 сутт T (GNED(2NI3) T `° T (QNED QR)... GNED 


Consequently the C, 's are determined and the function f(x) in real 


N 


time is 


a= > с 
=O 


VL 
on be Уш (58) 


Gauss-Laguerre Quadrature 
The Laguerre polynomials are the unique Gauss quadrature 
polynomials in the interval (0,9) with weighting function 
w(x) = Вт (59) 
Denote the Laguerre polynomials by LOO where 
Ы (х) = GDCD, (60) 
and the interval of expansion is (0,9). These polynomials are defined 


Ф 


by the Rodrigue's formula 


ох Eu (5 -ох) 
1 (ох) < е ак E е ; (61) 
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It follows from equation (14) that the orthogonality condi- 
tion is 
eo 
-QX a 
| e І. (ох) L, Cox) dx = 0 M+#N (62) 


О 


Consequently the resursion relations are 


еб айшы х) 1 (х) - КСО. (63) 


№1 
СӘКЕ ра оа Е (64) 


^ 


and 1 


The least squares polynomíal approximation to f(x) over (0,9) 


with respect to the weighting function 


w(x) = e”, (65) 
N 
is F(x) = > b L (ax) , (66) 
K=O 
where b. "DE | Ss f (x) L (ox) dx Ç (67) 
О 


Consequently the coefficients (or weights) become 


со L, (ax) 
-QX N 
A = | е c 1 (67) 
i (x x, ) Ly Cox, ) 
О 
Using the properties stated above, a similar procedure to that 


employed with the Legendre polynomials may be performed and an explicit 


inversion formula derived. 


TEST CASES 
A numerical technique has been developed based on shifted 
Legendre polynomials to invert the Laplace transform of a function. 
Let КЕ examine the accuracy obtained in utilizing this technique. 
Two test functions are considered. The first function is 
l 


u(t) = constant = 1 where the Laplace transform is v(s) = =e The 


; =t : 
second test function is u(t) = e | where the Laplace transform is 


T 


S +I ` The results obtained for these test cases are tabulated 


v(s) = 


29 


in Tables I and II. Notice that the accuracy decreases as N (number 
of roots) increases. However 8 - 9 figure accuracy is obtained in 
test cases I and 6 - 9 figure accuracy is obtained in test case II. 
These results indicate that the numerical inversion technique used 


is very accurate. 
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слее со пеши CXOQ C 


ОСОРО  ОСОООС ОСОС«СООӘОО 


Po Г САБЕ Т 
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290999399990 
SOLDI) 
49922299399) 


EOD 
“ОООО Ор 
a 29999999D 
2099 9999D 
коо д ОООО 
ые е е РЫУ 
ООЗЕ 


ee ee oD 
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91595999990 
25929 UID 
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52229999990 
27999999D 
So OD 
29299998992 
e299 9999) 
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.2999999997 
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„9999999992 
49999999907) 
999999999] 
.9999599991 
.999999999D 
92999999990 
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99999099990 
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00 
OO 
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OO 


CO 
QO 
QO 
OO 
OO 
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О.100000000р 
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0.1000000002 


0199959996560 
0220999996 в 
ООО ФЕ 
0.1000000032 
s. 090 oe” 
0.1000000С17 
099909 
0.999999944D 
0.1004000040 
98999985021) 
O.100000009D 
Ed D 
052209999620 
O.100000005D 


APPROXIMATE VALUE 


Ol 
00 
00 
00 


|= 


ШЕ) 


14 


ge CC 


i ee oe Cae 


er) ee CO OD) OO oe BO OG OE) OC. 


TABLE II 


EXACT VALUE 


10205681550 00 
.669990521D 00 
.330009478D 00 
,694318442D-0O1 


.97455 3956р 00 
20555920 09 
Ко 2225050 08 
712999999990 00 
.297077424р OO 
.1292344070 00 
» 2944004 36D-01 


.900953264D OO 
29253100D 0L 
203970475841 OO 
~710697697D ОО 
~974437169D OO 
.425;62830D OO 
22 720223020 00 
ШОО 65 90 
.074683166D-01 
.130407357D-01 


.993141904D 00 
.642174410 00 
ч ОПО ООО Т ПУ 
.043646452D 00 
го гогаз тор 00 
.659556184D OO 
коо ap 0 
e405 542525) OO 
БОЛЕЕ ОО 
s. 2 oo 00 
50555955 0D 
. 06 3993424D-01 
соого О 
.665009565D-02 
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АРРЗОХТМА Е VALUE 


079335601550 00 
0.669990521D OO 
0.330009478р 00 
0.6943184422-01 


0:974553056p 409 
0287076559208 
0.702922575D 00 
O. 500000000) СБ 
0.29707 74240) 090 
012923440 7 DM 
0.2544604 350-01 


0.95695 3264р 00 
0.93253165630 Ж 
0.9397047940 00 
271609769728 
0.574437169р 00 
0.4255628300 Oe 
0.2833023.2D O0 
0.1002952150 š 
0.0746631530200 
0.130467358D-01 


0.993141833р 00 
0.964217450D OO 
0.913600761D 00 
0.54364631Т7 00 
0.757024 350D age 
0.659556165D OO 
025540273990 099 
0.445972587D 00 
0.34044 3793D OO 
0.242315649D 00 
0.15635 3551004 
0.86 3993391D-O1 
0.357825607D-01 
0.685815140D-02 


III. THE INVERSE PROBLEM IN TRANSIENT HEAT CONDUCTION 


INTRODUCTION 
The problems presented in this section are representative of 
typical one, two, and three dimensional geometries. They were selected 
because of the existence of exact solutions, necessary for comparison 
of results. However, more complicated geometries can be accommodated 
by this method of solution, but the equations and analysis become 
quite complex. The FORTRAN IV programs used in the solution of these 


problems are listed in Appendix A. 


DIMENSIONLESS VARIABLES 
The solution to a problem in heat conduction can always be 
expressed in terms of dimensionless variables. This is usually done 
when making numerical calculations and is adopted here because of 
simplification in the differential equations involved and generality 
in application of the solutions obtained. 


The dimensionless variables used in the following problems are 


_ Kt Kt 
c EN De , 
d M ва 
се K = Cp(density) ' 
x r 
p Ер $ 
- uK uK 
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THE SEMI-INFINITE SOLID 
The method of numerical inversion of the Laplace transform 
is now applied to the solution of the semi-infinite solid. Since 
there exists no physical characteristic length, dimensional variables 
are used in the solution to this problem. 
The partial differential equation governing the temperature 


function is 


ou тои 

за NOx? ° e 
with the boundary conditions 

u(0,t) = Bt 3 (70) 
and 

u(x,0) 20 . (71) 


The solution of this direct problem is 


X 
р Kt 


Using this solution as the known temperature function at Xi (an 


(72) 





u(x,t) - 48ti^erfc 


interior location), we can solve the inverse problem to obtain u(0,t) 
at x » O (the surface of the solid). Here the exact solution is known 
as 

u(0,t) = ft. (73) 
t an aine the differential equation into the s plane we obtain the 


ordinary differential equation 


2 
57-%у-0 (74) 


where v is the Laplace transform of u(x,t) and is a function of x and 
s. The solution to equation (74) is of the form 
да ВИК x (75) 


v(x,s) 


Solving this equation at x = x, we find that 


nm vise FE Aj (76) 
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and therefore the solution of equation (74) becomes 


– 15 /К (x-x4) (77) 


v(x,s) - v(x,,s)e 


Since the temperature function at an interior location x, is given 


1 
in equation (72), we can solve equation (77) by taking the Laplace 
transform of equation (72) which then becomes 

v(x} >S) = Bre s/K ЕТ (78) 
Substituting equation (78) into equation (77) we have the solution 
in the s plane as 

v(x, 58) = (в, IS x) IE өсү) (79) 
At the surface of the solid x = 0, and the solution is 

v(0,s) s e „ЧЕЛ s). s/K zi (80) 
Taking the Laplace inverse of equation (80) we obtain the solution 
in the time plane at x = O (the surface of the solid). 

Àn illustration of the semi-infinite problem ís presented 
in Figure l. In obtaining the results plotted in Figure 2 it was 
found advantageous to use a convergence factor. The convergence 
factor e ^ was added to the function u(x,t) to improve the accuracy 
of solution. The function u(x,t) was multiplied by ole. transformed 
into the s plane and subsequently inverted back into the time plane. 
The final answer was then multiplied by e Note that the Laplace 


transform of this function is 


ife “u(x,t)} = Е(в+1). 


THE LONG CYLINDER 
The basic assumption involved in the solution of the "long" 
cylinder is that the cylinder is sufficiently long so that end effects 


can be neglected. In other words conduction is considered to be only 
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in the radial direction (Figure 1). The governing partial differential 


equation for the temperature function u(r,t) is 


pun Sow. 1900 
prm "= us зе) 2 ely 


subject to the boundary conditions 





u(R,t) = Bt , (82) 
u(r,0) O . (83) 


and 


In terms of dimensionless variables, the solution of the direct 


problem can be represented as the sum of an infinite series as 


2 a 2 
a 2. 1-p -ay т ___Чо(рд) 
u(p,tT) = Т 7 + пе е З J, ox) ) (84) 


where the Oy are defined by 
J (oy) =O . (85) 
Using this solution as the known temperature function at 04 
(an interior location), we can solve the inverse problem to obtain 
u(o,T) at p » 1 (the surface). Here the exact solution is known as 
Raat exe (86) 
Proceeding as before we find that the solution of the inverse 


problem in the s plane is 


I 
Du = v(p, »s) т. tp к (87) 





where 
EN Ер ла Jo(P19N) 
v(p4,s) с ти + 2 > Тр! Р Oy? J (On) J, (oy) ° (88) 


By numerically inverting equation (87), we obtain the solution of 
the inverse problem at the surface (o = 1) where the exact solution 
is given by equation (86). 

The results obtained for the "long" cylinder are plotted in 
Figures 5 and 6. It should be noted that as the distance from the 


surface is increased, the accuracy of solution decreases (Figure 5). 


This is reasonable since the temperature will be dampened more as the 


distance is increased. Let us consider a characteristic parameter 


l-r 
n- 


"E 


If a thermal penetration depth is taken as 

6 -2]|kt , 
then physical significance can be given to n. Now we can examine the 
accuracy of solution as n is varied. In Figure 5, surface temperature 
is plotted versus n for 6 = constant. It can be seen that as n is 
increased, the accuracy of solution decreases. Since the only variable 
in n is now r, it is obvious that as r increases, the accuracy of 
solution increases. By interpreting small n as "close" to the surface 


" from the surface, we can see the effects 


and large n as "far away 
of damping on the solution. 

A convergence factor was also used in this problem to improve 
the accuracy. As before, E was used as the convergence factor, and 
it can be seen from Figure 4 that the solution is improved. However 
the results obtained are not quite as good as the other problems 


considered. The inaccuracy introduced into this problem is due to 


the evaluation of an infinite series of Bessel functions. 


THE SPHERE 
The basic assumption involved in the solution of the sphere 
is axial symmetry. In other words the temperature function is not 
dependent on $ (Figure 1). The governing partial differential equa- 


tion for the temperature function u(r,t) is 


дъ 9?u , 2 дъ 
(2252 ы, (52) 
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subject to the boundary conditions 


ШЕ, Е) Eit (90) 


and 


u(r,0) 0 и (91) 


In terms of dimensionless variables, the solution of the 
direct problem can be represented as the sum of an infinite series 


as 
2 Ep a 2 cm =N? 7? 
u (09,1) » 1t - 2: - Tp í< N3 е ш T ein(Nmp) А (92) 








Using this solution as the known temperature function at 
01 (an interior location), we can solve the inverse problem to 
obtain u(p,T) at p = 1 (the surface). Here the exact solution is 
known as 
Mc T. (93) 
Proceeding as before we find that the solution of the inverse 


problem in the s plane is 
p; sinh (Js) 
К uet TER 


where 2 со N 
1 E 2 (-1) 1 i 
“(01.50 "3:57 7 тірі BE. Nš (amana o 


By numerically inverting equation (94), we obtain the solution of 


v(1,s) (94) 








the inverse problem at the surface (p = l) where the exact solution 
is given by equation (93). 

The results obtained for the sphere are plotted in Figure 7. 
Here as before a convergence factor си: was used to improve the 
solution. The results shown in Figure 7 are in very good agreement 


with the known exact solution. 
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INVERSE PROBLEMS IN TRANSIENT 
HEAT CONDUCTION 


SEMI- INFINITE SOLID 


ASSUMPTION ES 
CONDUCTION OCCURS 
ONLY IN X- DIRECTION 
7 
X 
X=0 X= X, 
SPHERE "LONG 
CYLINDER 
ASSUMPTION ASSUMPTION 
TEMPERATURE FUNCTION DOES END EFFECTS ARE 
NOT DEPENO ON d. THERE— NEGLECTED AND THUS 
FORE CONDUCTION OCCURS CONDUCTION OCCURS ONLY 
IN THE RADIAL DIRECTION. IN THE RADIAL DIRECTION. 


FIGURE | 
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u (0,T) (DIMENSIONLESS SURFACE TEMPERATURE. ) 


SURFACE TEMPERATURE vs. TIME 
FOR 
SEMI-INFINITE SOLID. 


å о 
о 
— EXACT 
^ NO CONVERGENCE 
FACTOR 


O CONVERGENCE FACTOR 


1 .e 25 .4 5 ‚6 0 


T (DIMENSIONLESS TIME. ) 


FIGURE 2. 
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и (О)/т) (DIMENSIONLESS SURFACE TEMPERATURE. ) 
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‚220 
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SURFACE TEMPERATURE vs. TIME 


.04 


FOR 
"LONG" CYLINDER 


08 12 16 
T (DIMENSIONLESS TIME) 


FIGURE 3. 
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.20 


— EXACT 
o P 

^ P -.5 
D P =1.0 


‚224 


‚28 


u (0,t) (DIMENSIONLESS SURFACE TEMPERATURE.) 


‚32 


‚28 


‚24 


‚20 


. ОВ 


SURFACE TEMPERATURE vs. TIME 
FOR 


"LONG" CYLINDER 
( CONVERGENCE FACTOR ) 


—— EXACT 


O APPROXIMATION 
WITH f, =! 


A APPROXIMATION 
x р, = .8 


.04 .08 12 16 ‚20 ‚24 ‚28 


t (DIMENSIONLESS TIME) 
FIGURE 4. 
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u (O,T) (DIMENSIONLESS SURFACE TEMPERATURE.) 


‚32 


.28 


‚24 


‚20 


. 08 


.04 


SURFACE TEMPERATURE vs. TIME 


FOR 
SPHERE 
^ 
© 
— EXACT 
^ APPROXIMATION 
О APPROXIMATION WITH 
CONVERGENCE FACTOR 
.04 .08 42 16 .20 ‚24 ‚28 


t ( DIMENISONLESS TIME ) 
FIGURE 6. 
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DISCUSSION OF RESULTS 

In the determination of surface temperature by the method of 
numerical inversion of the Laplace transform, it was found necessary 
to utilize a convergence factor to improve accuracy. The function 
u(x,t) was multiplied by Са transformed into the s plane, and 
numerically inverted into the time plane. Then the function obtained 
was multiplied by еб to obtain the desired surface temperature function 
u(0,t). This method caused the approximation to converge to the known 
exact solution. In the semi-infinite solid and the sphere, the 
results obtained from the use of the convergence factor were very 
good. In the "long" cylinder the results obtained were not quite as 
good. The inaccuracy introduced in this problem is due to the evalua- 
tion of an infinite series of Bessel functions. 

The damping effects were considered for the case of the "long" 
cylinder. It was found that as the distance from the surface was 
increased, the accuracy of the solution decreased. This was a 
reasonable result, since it was logical to expect the solution to 
lag as the distance from the surface was increased. 

The use of time scaling (mentioned in section II) was only 
incorporated in the semi-infinite solid FORTRAN IV program, but it 
can be included in the other problems. Although the initial tempera- 
ture distribution was taken as zero, any reasonable initial tempera- 
ture distribution can be handled with this method. Such a change 
would simply result in more terms to be inverted. Also it should be 
noted that although the solution was obtained for the surface tempera- 


ture, heat flux at the surface can also be determined. 
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CONCLUSIONS 

Three geometries were considered in the solution to the 
inverse problem in transient heat conduction. These geometries, 
the semi-infinite solid, "long" cylinder, and sphere, are typical 
of one, two, and three dimensional problems respectively (it should 
be noted here that an analogous approach can be made in the solution 
to the finite plane bounded by two parallel planes). Solutions were 
obtained which agreed very well to the known exact solution. 

Several parameters were investigated with respect to the 
effects on the accuracy of solution. It was found in all cases that 
utilizing a convergence factor caused the solution to better approxi- 
mate the known exact solution. The damping effect was considered in 
the "long" cylinder problem. As expected it was found that as 01 
increased (distance from the surface decreased) the accuracy of 
solution increased. 

Therefore the solutions obtained indicate that the method of 
numerically inverting the Laplace transform provides an accurate 


approximation to the known exact solutions. 
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ту. SUMMARY 

The numerical technique involved in numerically inverting the 
Laplace transform was methodically developed. Elementary theorems 
of the Laplace transform were introduced, and the instability of the 
inverse of the Laplace transform was briefly discussed. Gaussian 
quadrature methods were subsequently introduced and developed. The 
shifted Legendre polynomials were then examined for several reasons. 
It was found that the coefficients of the shifted polynomials were 
integers, and the region of orthogonality (0,1) was found to be most 
convenient in the following development. Then, an explicit inverse 
matrix was found which would be used later to solve the inverse prob- 
lem in transient heat conduction. Two useful numerical properties 
were discussed. Interlacing of the roots (and therefore times) was 
found useful in interpolation of the results, while change of the 
time scale was found useful in ER the range of times. Two 
alternate methods of solution (Papoulis and Gauss-Laguerre quadrature) 
that could be used to obtain an explicit inverse matrix were subse- 
quently introduced. Finally, the explicit inverse matrix was tested 
for accuracy. The two test cases employed were for the functions 
u(t) = 1, and u(t) = ai The accuracy was found to be very good for 
both cases. 

Three geometries were considered in the solution of the 
inverse problem in transient heat conduction. These problems, the 
semi-infinite solid, "long" cylinder and sphere, are typical of one, 
two, and three dimensional problems. They were chosen because of the 


existence of exact solutions. It should be noted that the results 


obtained were for discrete points and not for a continuous spectra of 
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information. Consequently this method is applicable to problems 
where discrete information is known (i.e., data from experiment). 
Several parameters were examined with respect to influence 
on the accuracy of solution. A convergence factor eda was found 
to improve the accuracy of solution in all cases considered. It 
was found that as the distance from the surface was increased, 
the accuracy of solution decreased. This was considered reasonable 
since damping effects should be more noticeable at greater distances 
from the surface. 
Consequently the solution to the inverse problem in transient 
heat conduction by numerical inversion of the Laplace transform was 
found to yield accurate results for semi-infinite solid, "long" 


cylinder, and sphere. 


44 


OF 


E. 


12. 


D: 


BIBLIOGRAPHY 


Beck, J. V. "The Inverse Problem in Transient Heat Conduction," 
Journal of Applied Mechanics (June, 1965), pp. 472-473. 


Bellman, Richard, Rebert E. Kalaba, and JoAnn Lockett. Numerical 
Inversion of the Laplace Transform. New York: American 
Elsevier Pub lishing Company, Inc., 1966. 


Bellman, Richard, Robert E. Kalaba, and Marcia C. Prestrud. 


Invariant Imbedding and and Time-Dependent Transport Processes. 
New York: American Elsevier Publishing Company, Inc., 1965. 


Burggraf, О. R. ‘An Exact Solution of the Inverse Problem in 
Heat Conduction--Theory and Applications," Journal of Heat 
Transfer (August, 1964), pp. 373-382. 


Carslaw, H. S. and J. C. Jaeger. Conduction of Heat in Solids. 


2nd ed. London: Oxford University Press, 1959. 


Choquet-Bruhat, Yvonne. Problems and Solutions in Mathematical 


Physics. Translated by C. Peltzer. San Francisco: Holden- 
Пазл Inc... 1967. 


Cost, Thomas L. "Approximate Laplace Transform Inversions in 
Viscoelastic Stress Analysis," AIAA Journal (December, 1964), 
pp. 2157-2166. 


Cost, Thomas L. and John J. Brisbane. Quarterly Progress Report 
on Engineering Research (U). Report No. P-63-13. Huntsville, 
Alabama: Rohm & Haas Company, Redstone Arsenal Research 
Division, 1963. 


Costello, Frederick A. and George L. Shrenk. "Numerical Solution 
to the Heat Transfer Equations with Combined Conduction and 
Radiation," Journal of Computational Physics, Vol. 1l, No. 4 
(June, 1967), pp. 541-543. 


Crank, J. The Mathematics of Diffusion. London: Oxford 
University Press, 1956. 


Deverall, Г. I., and R. S. Channapragada. "A New Integral 
Equation for Heat Flux in Inverse Heat Conduction," Journal 
of Heat Transfer (August, 1966), pp. 327-328. 


Frank, Irving. "Ап Application of Least Squares Method to the 
Solution of the Inverse Problem of Heat Conduction," Journal 
of Heat Transfer (November, 1963), pp. 378-379. 


Jahnke-Emde-Losch. Tables of Higher Functions. Revised by 


Dr. Friedrich Losch. New York: McGraw-Hill Book Company, 
Inc., 1960. 


45 


14. 


lie 


16. 


T77 


To: 


О 


20 


Klerer, Melvin and Granino A. Korn (eds.). Digital Computer 
User's Handbook. New York: McGraw-Hill Book Company, 
Lnc оО 


Krylov, Vladimir Ivanovich. Approximate Calculation of Integrals. 
Translated by Arthur H. Stroud. New York: Macmillan Company, 
1962. 


Lanczos, Cornelius. Applied Analysis. Englewood Cliffs, 
New Jersey: Prentice-Hall, Inc., 1956. 


Saaty, Thomas L. Modern Nonlinear Equations. New York: 
McGraw-Hill Book Company, Inc., 1967. 


Sparrow, E. M., A. Haji-Sheikh, and T. S. Lundgren. "The Inverse 
Problem in Transient Heat Conduction," Journal of Applied 
Mechanics (September, 1964), pp. 369-375. 


Stolz, G., Jr. "Numerical Solutions to an Inverse Problem of 
Heat Conduction for Simple Shapes," Journal of Heat 
Transfer (February, 1960), pp. 20-26. 


Stroud, A. H. and Don Secrest. Gaussian Quadrature Formulas. 
Englewood Cliffs, New Jersey: Prentice-Hall, Inc., 1966. 


46 


APPENDIX A 


FORIRAN IV PROGRAMS 


]. The Semi-infinite Solid 
The following program is listed: 
SEMI-INFINITE PLANE SOLUTION USING DIMENSIONLESS 


VARIABLES IN ONE-DIMENSIONAL SYSTEM. 


The following subroutines are needed: 


INPUT Subroutine 


EXPINV Subroutine 
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FORTRAN IV PROGRAM SYMBOLS 


MAIN PROGRAMS 


Т1 
ТАП 
RO 
ROL 
APPR 
XK 
EXACT 


SCALE 


ALFA 


ARQ 


BESSO 
ВЕ551 
BISS 

BISSI 


ADA 


Explicit inverse matrix 
Roots of the shifted Legendre polynomials (х.) 

х 
Weights of the shifted Legendre polynomials (А) 


Values of time (t; = -log х.) 


Non-dimensional time (T = БЕ - Zr 
Distance (p) 
Interior location (ру) 


Approximate value of v(x,s) 


Reciprocal of K 


Known exact solution 

Scale factor (a) 

Value of ст: UT) 

Roots of the equation J (G) - 0 
Argument of the function J Сок?) 
Argument of the function I (ds 04) 
Argument of the function I (ds ) 


Resultant value of J Cop) 


of J4 Co) 
of I (Qs ) 
of I (ds 04) 


Resultant value 
Resultant value 
Resultant value 
Value of n 


Number of roots 
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EXPINV SUBROUTINE 


C Legendre coefficients 

* 
B Coefficients of the polynomials Р(х) /(х-х.) 
Е Legendre derivatives 
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2. The “Lonzi ylinder 
The following program is listed: 


TWO-DIMENSIONAL "LONG' CYLINDER SOLUTION. 


The following subroutines are needed: 
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3. The Sphere 
The following program is listed: 


THREE-DIMENSIONAL SPHERE SOLUTION. 


The following subroutines are needed: 
INPUT subroutine 


EXPINV subroutine 


DS 


SOLUTIQN ПЕ THE 
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4. Subroutines 


The following subroutines are listed: 


INPUT subroutine 
EXPINV subroutine 
* 
BESI subroutine 
* 
BESJ subroutine 


ES 
"Subroutine BESJ" and "subroutine BESI," System/360 


ee eS 


Manual. White Plains, New York: International Business Machines 
Corporation, 1967, pp. 156, 157, 159; 
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SUBROUTINE RESJ әд 


PURPOSF 
| COMPUTE THE J BESSFL FUNCTION FOR A GIVEN ARGUMENT AND ORDE 


.- 


USAGE 
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AND M. ABRAMOWITZ,aGENERATION OF BESSEL FUNCTIONS ON HIG | 
SPEED СОМРИТЕВ 52, М.Т.А.С. » 1111957, РР. 255-257 


геегееегееееееееееееееееееееееееегеееееееееееееееегееееегеегегееееееее 


SUBROUTINE RESJIX Ns BJ D, IFR) 
EE DU n e NE 


Gu yD CCDC IC GC ICO IO ICO IT IOI OIC IO IMI INN 


-- 
2 
WZ LE 
> lo 
е 
UJ 
< 
a 
UJ 
є— 


$ 29 9 22 | n MEME қ > 
x 
x 


| 
| wv 
12 2 


Len E s E SSS SU r u 1. = CNN O a и 


OJ UJ OJ WNW aJ 
j co^ Nr 

Z ZM 
aA OWN NK 


I 
* 236^! A04! Xm—Z2x-c-—O 


Ics ee a) es T) eat ec ZU ec D 
"Ur-mrma^nT-AO-nnm тїгєт TNL 
ох! 2-4--11 гп 4 20 2 420 
"mz Си 200 Чо Н >< C Il Ze 


, 
Eel STARTINGSVADUESORSM—————————— M B 1 


ТЕ (Х-5. ) БС „60, 60 | 
50 MAzIiX*6, я | 
ато а ee | 
60 MA-l. 6X *60, /X | 

70 MRZzN*IFIX(X)/4-42 
м7 ЕВО= МАХО (МА, МВ) | 


-— ` s + . Cor + oe a “— m pa qr  — —  — eee e ER — ә - . --- — 2 


m 
С SET UPPER LIMIT OF M 
MMAX=NTEST __ 


Р 


| 
AAO 








6Q 


—F n MMMÀ  — = — —— —_ <= - - ` _ - — =s 


150 DO 19^ M 


M7FRO+MMAX 93 


а RR — - = <+ - 


SET F(M) ,F(M-1) 


CuO 


2) 2) ТЕГО" lee 


0 
/ 


LL LL — — CS ae сички 72%7 


22545 


— NIA 
=—t rot od 


FM] 
м1 = ВМК 


---- ee .. 


Ge ND нас па аа 


ALPHA#RMK*S 


2, ЖЕМ| /Х-ЕМ 


+ 

саво 
Н 

В 


О До т 


Мк 
F 


ALPHA*BMK 
J/ALPHA 


MK 
As 


IF(DARS(BJ-BPREV)-DABS(D*BJ))200,205,190 


190 BPREV=BJ 


TER=3 
200 RETURN 


END 


- — —— —— —xÑII:- 


— — — ee ate ee a ee ты = 


— — —- a rn ee ee ee en we m os - 


` 
——.— = = am 


-—— — ee 


— — а 


-m — — — — —— — ыыт * 


6 


esl 


APPENDIX B 


ROOTS OF THE SHIFTED LEGENDRE POLYNOMIALS 
AND CORRESPONDING WEIGHTS 


* 
The roots, Xi» of the shifted Legendre polynomial PN (x) are 


computed from the roots of Legendre polynomial Р (г) using the formula 





The coefficients (or weights) of the shifted Legendre 


* 
polynomial Ры (x) in the interval (0,1) are then 


* 1 
Ау ри , 


where А, 5 are the weights of the Legendre polynomial Pr) in the 


interval (-1,1) and were obtained from Stroud and Secrest [20]. 
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ROOTS AND WEIGHTS 


ROOTS 


0.788675134594812600 00 
0.211324865405186900 00 


ощсстезовзаоготаперо 00 
0.5 00000000000000000 00 
O 112760 653710 5300 00 


0.93056815579702600р 00 
0.60999052179242750D ОО 
0.33000947820757180р OO 
0.694318442029738000-01 


O" 0522965) I20D 00 
0.769234655052841300 00 
0.5 00000000000000000 00 
0.230765 34494715830р 00 
0.46910077030668140D-01 


0.966234757101575700 00 
0.830604693233131900 00 
0.61930959304159840D 00 
0. 38069040695840140D 00 
0.169395306766867400 00 
0.337652428984240400-01 


O .97455395617137920D ОО 
r O7076559279969690D OO 
0.70292257568869830D 00 
0.500000000000000000 00 
0.29707742431130150D 00 


0.12923440720030290D OO . 


0.254460438286207700-01 
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NEIGH?S 


0.5 00000000000000000 00 
0.5 00000000000000000 00 


0.2777777777771777700 00 
0.444444444444444200 00 
0.277777777777777700 00 


0.17392742256872660D 00 
0.32607257743127270р 00 
OS 3260725714 M72 р оо 
0.173927422566726600 00 


0.118463442528094400 00 
0.239314335249682900 00 
0.284444444444444300 00 
0.239314335249682900 00 
0.118463442528094400 00 


0.85662246189584 8000-01 
0.18038078652406920р 00 
0.233956967286345 30D 00 
0.23395696'7286 345 30D 00 
0.18038078652406920D 0C 
0.856622461895 84800р-01 


0.647424830844345000-01 
0.13985269574463800р 00 
ша еще 025295 5591 D 00 
0.208979591836734300 00 
0.19091502525255910D OO 
0.13985269574463800D OO 
0.647424830844 34500D-01 


ROOTS AND WEIGHTS (continued) 


ROOTS 


0.98014492824876800D 00 
О.89833323870681 3001 00 
0.76276620495816400р 00 
0.59171732124782460р 00 
0.40828267875217520D 00 
0.297233 719509 «*195550D 00 
0.10166676129318050р 00 
0.19855071751231980D-01 


0.98405011975 3812700 00 
0.91801555 3663317400 00 
0.50668571635029520р 00 
0.662126711701904 300 00 
0.500000000000000000 00 
0.337873288298095500 00 
0.193314283649704600 00 
0.819844463366819000-01 
0.15919880246187050D-01 


0.98695 3264258585401 00 
4.93253168334449190D OO 
0.53970478414951200D OO 
0.716697697064623500 00 
0.57443716949051550D 00 
0.42556283050915410р 00 
0.283302302935 376300 00 
0.160295215850487500 OO 
0.674683166255077000-01 
0.130467357414141900-01 


0.989114329073028200 00 
0.94353129958404730D 00 
0.86507600278702430D 00 
0.75954806460340560D 00 
0.63477157797617220D 00 
O.50000000000000000D 00 
0.365228422023827300 00 
0.24045193539659390D 00 
ова троен е9тоео0 00 
0.56468700115952460D-01 
0.108856709269715200-01 
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10 


WETO TED 


0.506142681451880600-01 
0.111190517226686900 00 
0.15685332293894350D OO 
0.18134189168918070D OO 
0.151341891689180700 00 
O.15685332293894350D 00 
0.11119051722668690р 00 
0.50614268145188080р-01 


0.40637194180787150D-01 
0.90324080347428300D-01 
0.13030534820146730D 00 
0.15617353852000140D 00 
0.16511967750062960D 00 
0.15617353852000140D 00 
0.13030534820146730D 00 
0.90324080347428300D-01 
0.406371941807871520-01 


0. 33335672154344030D-01 
0.74725674575290000р-01 
0.109543181257990800 00 
0.134633359654997800 00 
O.14776211235737600D Os 
0.147762112357375000 8 
0.134633359654997800 00 
0.109543161257990800 00 
0.74725674575290000р-01 
0.33335672154344030р-01 


0.27834283558086790р-01 
0.627901847324519000-01 
0.931451054638668000-01 
0.116596822295995100 00 
0.131402272255 1229001068 
0.13646254338894990р 00 
0О.131402272255 122900898 
0.11659688229599510D OO 
0.93145105463866800р-01 
0.627901847324519000-01 
0.276342635530867900-01 
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ROOTS AND WEIGHTS (continued) 


ROOTS 


.99078031712335930D 00 
.952055628155237000 00 
.88495133709715200D 00 
.79365897714330040р 00 
.68391574949908970D 00 
.56261670425573420р 00 
.43738329574426560р 00 
.31608425050090960D 00 
.206341022056091100 00 
.11504866290264750D 00 
.47941371814762010D-01 
.921968287664042501-02 


.99209152735929360р 00 
.95879919961148360D OO 
.900789045360665450D OO 
.02117466972016980D OO 
2A246329 11022310D 00 
Molo 22915797 (56 730n. 00 
.50000000000000000р 00 
.384770842022432300 ОО 
"7575562 4451776401] ОО 
» 178825 33027952970D 00 
. 99210954633 345100р-01 
.41200800388511100р-01 
.79084726407060300р-02 


.99314190434840570р 00 
.96421744183178630р 00 
.91360065753488210D 00 
.84364645240584240D 00 
.757624318179076Т700 00 
.659556184463944600 00 
.55402747435367150р 00 
.44597252564632830р 00 
.34044381553605520D 00 
.24231568182092280D 00 
.15635354759415710р 00 
.86399342465117200D-01 
.35782558168213250D-01 
.68580956515938870D-02 
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12 


13 


eee OO OOOO 


о аа оо ооо оо 


оосо ОСЕ СЕСЕ ОЈО 


WEIGHTS 


.23587668193255860D-01 
.534696629976591800-01 
.80039164271672700D-01 
210159371 2 ЕЕЕ ОО 
.11674626826917710D 00 
.12457352290670110D 00 
.124573522906701100 00 
.11674626826917710D 00 
.10158371336153270D 00 
.80039164271672700D-01 
5 346960299765 9150D-01 
-2 3587668193255 860D-01 


» 2024 2002 382657900D-01 
.460607499188641707-01 
-694 3675510959 3600D-01 
-59072990 38097 2600D-01 
.10390802376844400D 00 
.11314159013144530р 00 
„1162 ОЕ p oO 
„11.31415 9048542 поцатиоо 
.10390802376844400р 00 
. 59072990380972600р-01 
.69436755109893600р-01 
.46060749918864170D-01 
.20242002382657900р-01 


.17559730165875880р-01 
4007904 3579880060D-01 
60759285 343951540D-01 
7860158 3579096500D-01 
927691987 38968500D-01 
.10259923186064770D 00 
.10763192673157860D 00 
.10763192673157860D 00 
.10259923186064770D 00 
9276 91987 38968500D-01 
. 786015 83579096500D-01 
60759285 343951540D-01 
4007904 3579880060D-01 
175597 30165875880D-01 


о. 
O. 
O 
.86220580563008460р 00 
.785486005630426920р 00 
-69707567353875130D 00 
.60059704699871720р 00 
.50000000000000000D 00 
. 39940295 300128260D 00 
.302924326461219000 00 
„214513913695730600 00 
ТЕ ЦИ ЗТОФИЯТОТ OO 
.75 896708294786100р-01 
.31363303799647090D-01 
.60037409897574110р-02 
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ROOTS AND WEIGHTS (continued) 


ROUIS 


99399625901024220D 00 
96863669620035250D 00 
ОСАО роза арр OO 
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15 
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WEIGHTS 


.15376620998056590D-01 
© 35183023744054020D-01 
.5 357961023 35 85930р-01 
.69785 3 38963076900р-01 
.83134602908496700р-01 
- 930805000077 80800D-01 
.99215742663555600р-01 
.10125912096278020р 00 
9921574266 3555600р-01 
.93060500007780800р-01 
.531346029084967000-01 
.69782 33896 3076900р-01 
.535196102335959300-01 
.351930237440540200-01 
.153766209950585907-.1 


APPENDIX C 
Legendre Coefficients 


This table gives the coefficients of the first fifteen 
* 
shifted Legendre polynomials, Ры (x), in order of lowest to highest 


power such that 


со 
* K 
Pa (x) = > ES Е 


к<о 
where P. (x) = Р,(1-2х). Тһе recurrence relation for these 
coefficients is 
(КОС < (28-30, 1, - сш - Снори и : 
All of these coefficients are integers but due to machine conversion 
from hexadecimal to decimal, noninteger numbers may be printed in 
this table. In other words 
-0.13999999999999990D 03 represents the integer -14. The coef- 
ficients are printed in order of lowest to highest power of x. For 


example 


* 2 
Р, (x) = 1-6x-6x'. 
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-0. 


-0. 
O. 
-0 
О 


LEGENDRE COEFFICIENTS 


N= 1 


0999999920980 900) 
.200000000000000000 


N = 2 


«ESO IM OD 
„2 ОООО OD 
2999992 89999880 


N = 3 


soo eee 99999900) 
-119999999999999980 
299989999999993960р 
. 200000000000020000 


N = 4 


„заеззазаввезазчвор 
. 20000000000000000D 
-99999999720 J9990D 
«13999990399790 
.70000000000000000D 


N= 5 


ss MOD 
saco 
.21000000000000000р 
.55999999999999980р 
5=6с53209995995399508 
epee Ор 


N = 6 


99999999893 9900D 


fe 9999999999 70D 


Oma 9999999999990) 


5072 899999999900 
31499999999999990D 


2217125292099 99999990D 
seme 79999999990) 


OO 
O1 


OO 
01 


OO 
O2 
02 
O2 


OO 


02 
03 
02 


OO 
02 
03 
03 
03 
03 


00 
02 
03 
04 
04 
04 
03 
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N = 7 


Qa89099999999999909D 
-019559999999598080090 
0.756000000000000000 
-O 4419999999 99S ши 
0.11542999999999990D 
-0.16631999999999990D 
Ogl201l99 gode 
2963431999 9099959 19) 


N-8 
0.99999999999999900D 
-0.720000000000000000 
0.12599999999999990D 
-0.923:9999999999990D 
0.34649999999999980D 
-0.72071999999999980D 
0.84083999999999980D 
-0.51479999999999990D 
0.12869999999999980D 


N = 9 


0.9999 JE E ROCCO] 
-O.O9 5 j D 
0.198000000000000000D 
-0.154799999999999980D 
0 < 900r 99999600 
=O0:2 5109993995 55 ЕЕ 
ОСЕ ОА 0099995993680 
01-4 1019559999999 2859090 
0.210789999999999980D 
-0.406199999999999970D 


N = 10 


O. 9599999999 eu 
-O.1050999999999 5 imn 
O: 6999999 s 
-02945 199999999 IM 
0.210209999999999960D 
-0.126797299999 DIM 
О -ш6526/7 9999998989980 
Ос 5 5 OD 
0.196910999999999900D 
-0.923779999999999980р 
О а j OD 


LEGENDRE 


N - 11 


0.99999999999999900D 
-0.131999999999999907 
0.42599999999999950р 
-0.000599999999999800 
0.45044999999999980D 
-0.201601599999999501 
0.57177119559999990р 
-0.105.71919999999990р 
0.124710299999559500 
-0.92377999999999980р 
0.38798759999999970D 
-0.70543199999999980D 


N= 12 


0.99999999999999900D 
-0.155999999999999900 
0.60059999999999980D 
-0.1С009999999999980р 
0.90089999999999980D 
-0.490089599999999901 
0.17153135999999990D 
-0.39907295999999580D 
0.62355149995999980 
-0.64664599999599980D 
0.426786 35999999990D 
-0.16224935999999990D 
0.270415600000000001 


Мя 13 


oE 99 999999900D 
со ша е зсачвзчззор 
ОХО 8892 29 39999999У9Ср 
2157501509999 9999960D 
PE ROS EU OUgg99oOD 


-0.110270159999999500 


0.46558511999999980D 
-0.13302431999999 90D 
0.26189162999999950D 
-0.35565529999999980D 
0.327202575999993900 
-0.194699231999999800 
0.67603899999999980р 
-0.10400599999999980D 


COEFFICIENTS (continued) 


OO 
03 
04 
05 
06 


O7 
ere 


ОТ 
07 
06 


OO 
03 
04 
06 
06 
07 
09 
09 
08 
08 
08 
08 
07 


00 
03 
04 
06 
ОТ 
08 
09 
e 
09 
09 
o9 
09 
09 
08 
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N= 14 


.99999999999999900D 
.21000000000000000D 
.10920000000000000D 
2415200000000 C 00C D 
- 306 30600000000C00D 
.23279255999999980р 
.11639627999999980р 
. 39907295999999980D 
.900269309999999900 
.163601437999999907 
. 196 32172559999480D 
.16224935999999990D 
.87585069999999990р 
.28081619999999950р 
.40116599999999990D 


N= 15 


29222929999 В ЕО 
3997999959 UID 
.14279999999999980р 
«3127909900 E 
„2290739992 nD 
.46558511999999960D 
42112013 Ор 
„10974506990 р 
» 31551705899999980D 
.654405751999999700 
.991609627999999600 
.10546208399999980р 
. 19096562999999960D 
e 3 19142679 0 9970D 
.11633813999999982р 
vues БЕ 


APPENDIX D 


Negatives of Logarithms of Roots of 
Shifted Legendre Polynomials 


The negative logarithms computed in these tables are 


t, 
1 


Р (х.) 


-log x, , 
1 
where 


D. 3 i h. 12. 0 is 
N 


ii 
(9 
+ 
> 
мә 
хә 
H= 
Ul 


It can be readily observed from these tables that the times 
or t values interlace for different N, making possible the use of 


more t values for extrapolation purposes. 
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соооооО Caro 2 © GEO Coe 


Ore CO 


О 
О 
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NEGATIVES OF LOGARITHNS OF ROOTS OF SHIFTED 
LEGENDRE POLYNOMIALS 


N=2 


.237400786151601920р 00 
.155435060907502999D “Ol 


N = 3 


.11957401204924280D OO 


.69314718055994510р 00 
.21830110809448010р 01 


. N = 4 


.719599591587634000-01 


.40049171327565440D OO 


.11086339030929460D O1 


1266749996665219690D 01 


М = 5 
.48046021994112870р-01 


зер зешедллобберр ОО 
Сакс 85 1545100 00 


.146635 39074975080р 01 
„305952276454629200 01 


N = 6 
. 3434645452025 2990D-01 
.18560129738533360D O 
.47914995102148920р 00 


.965776801433932520D OO 


SUE D20202017DAT ЗОР 01 
шро во дер стогов 01 


М-7 


.25775393510447510р-01 
Meroe 46259 94 3310D 00 
. 352508527 39216200D 00 
.69314718055994510р 00 
.12137624562357850р 01 
E204 Cl1274134878730D OI 
.36711949968461940р 01 


-LOG(x) 
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CO C O OO CGO O Oe Qe OOO 


@ © CO OC CO Oe @ 


O O C Фф д» ею еы со ар ер 


М = 8 


.200548 322834406 30D-O1 
.10721418964279870р 00 


.270803710172606800 
.52472525602874090р 
.89579550440199850D 
14 3870914712588 30D 
.22860548602906120D 


OO 
OO 
OO 
91 
91 


. 39192958006440780р 01 


М = 9 


.160479627304791900-01 
.8554094551698139000-01 
«21452113346 20C T ODTOD 
.41229833395649940D 00 
.69314718055994510В 00 


„1085084 3404 7872601 
.164 34 3800181081 30D 
250122571995 39] 590 
.414018662081862000 


N = 10 


.1313259190609151910В- 
.698521512644286000- 
.174704896357960300 
.33310111 ЕШ OD 
.55436455322540930р 
. 354 34267909849090р 
СОБЕ ТОА ОША аро 
ОД ре а 4/ uM On 
.269609717417964980D 
.43392173112598770р 


N = 11 


.10945353360154720D- 
.58125740495859540D-01 
.14493791141644620D 
.21503167441491130D 
.45449006452666760р 
.69314718055994510р 
.10072323073755050р 
.14252350622303700р 
200304 36426860710D 
.28740687745821190р 
.45203079484393840р 


Ol 
01 
on 
Ol 


01 
01 
ОО 
ОО 
OO 
OO 
Ol 
Зи 
О1 
Ol 


OU 


OO 
OO 
OO 
00 
01 
01 
Ol 
Ol 
Ol 


NEGATIVSS-.OF LOGARITHMS..O# ROOTS OF SHIFTsD 
LEGENDRE POLYNOMIALS (continued) 


N = 12 


0.92624472047904000D- 


O2 


0.491266618588122200-01 


0.122922929/ се а ау) 
0.23110140981597780р 
0.37992054219297830р 
0.57515669219094 1000 
0.026945 36135 380250D 
0.115174648542161401 
О По So2250201 3G eon 
0.21624000844470620D 
O. 30377764 351420280D 
0.46864146 3716 34650D 


N = 13 


0.79399104703031280D- 


00 
00 
00 
00 
OO 
Ol 
Ol 
Ol 
Oil 
Ol 


O2 


“= 4 O O 087655600=01 


0.10448418267570260D 
0.19701943975857120D 
0.32262364680743880D 
0.48576046598152780D 
0.69314718055994510D 
0.95510733743158590D 
0.12882474768486800D 
0.17213457580450050D 
0.23105068410166020D 
0. 31892975959097460D 
0.48398206080462550D 


OO 
OO 
OO 
00 
00 
00 
01 
Ol 
ol 
01 
Ol 
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N= 14 


0.69817204656915940р-02 
0.364384477407208900-01 
0.90361720362083900D-01 


S CGO OOOO O O O 


ао Ооо о СО оО 


.17002176740104690р 
.27756763865621490р 
.41618811793798540D 
.59054100076837180р 
.80749793054938170р 
.1077505172775 3410р 
.14172663523405600р 
.18556355052547450р 
.244877521355851800 
.33302947063021520р 
.49823254775913270р 


М = 15 


.602183590382997302- 
. 31865653915715830р-01 
.18931426028461700D- 
. 1482577 3409602940D 
.241452534614226900 
. 36086130375801570D 
.50983104019704670D 
.69314718055994510D 
.91778446443415 1100 
.119427225231764T70D 
.15 39 3806770033880D 
. 198201625977 33620D 
2578381964496 1590D 
. 34621167 382948760D 
.51153725057554800р 


00 
00 
00 
00 
00 
01 
Ol 
OT 
o 
Qu 
Ol 


02 


01 
00 
00 
00 
00 
00 
OO 
Ol 
Ol 
Ol 
Ol 
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Ov 


APPENDIX E 


Legendre Derivatives 
The derivatives of the shifted Legendre polynomials tabulated 
in this table are evaluated at the roots and listed in order of 


smallest to largest root. For example, 


ж! 
P. (0.5000) = -3.75. 


Recalling that the recurrence relation for the derivative of 


the Legendre polynomial is 


CED P COM 


N 


and the relationships between shifted Legendre polynomials and Legendre 


(х) = м ХР (x)-7Py 1 


polynomials are 


* 
Py ex) = Р „(1-2х) к 
апа 


NN 2Р (1-2 
we can determine the derivative of the shifted polynomial from the 


formula 


d * ж * 
-2x(x-l)4- Py (x) = N[(1-2x)P,, (x)-P& 160] : 


T3 


-0. 
90 р 
.99999999999999990D 


-0 


-0 


-O. 
.485177807260092500 
-0. 
.485177807260083801 
-0. 


-0 


=" 


-0 


-0 
-0 


LEGENDRE DERIVATIVES 


N23 
5999999999009 CRISI» 


N = 4 


.94 353 7565  370900D 
.37243057339114480D 
.37243057339114550D 
-0. 


94332756589371720D 


М = 5 
13740666961459280р 


SA oe 299990950) 
13740666961489610D 
N = 6 


CESSIT S504 156095 
Hoe (Ono Leo I OPEN 
.42578674578469820р 
-0. 
.6277070185998644000 
.1891598687 3542170р 


42578674578469820р 


N = 7 


.24956992249836780р 
. (3 /12-2159-4960360D 
=. 
.43149999999999950D 
-0. 
-197122139249200500 
-0. 


50083119828975910D 
50083119828975930D 
24956992249837180D 


N= 8 


.31862704825655890D 
-0. 
59356616715 355440D 
47776 387675098840D 
47776 387675098840D 
Пи оо ЗОП 
.99233405513433910D 
E 


992334055134 32490D 


318627048256566 30D 


DJ 
01 
91 


02 


02 
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м+9 


-0.39632641191048770D 
0.121285105985909201 
-0.701511394936625700 
0.534993951976716500 
-0.49:187499999999801 
0.534993851876717201 
-0.7015113984936626800 
0.12128510598591260D 
-0. 396 32641191052960р 


N = 10 


0.482665 3972945 8070D 
-0 14584224552 352420D 
0.8235 3816679841920D 
~0.604826 20784 350310D 
0.52615698870016580D 
~0.52615698870016690D 
0.60482620784 3503600 
-0.82353816679842400D 
0.14584224552352820D 
-0.482665397294655700 


№ = 11 


-0.57764248669799700р 
0.17269087349651640D 
-0.95906545604470890D 
0.685274800377843500 
-0.57293788755183090р 
0.54140624999999990р 
-0.57293788755183200р 
0.68527480037704640D 
-0.959065456044717600 
0.17289087349652200D 
-0.577642486698071501 


02 
O2 
01 
en 
DI 
OL 
Ol 


02 


02 
02 
01 
Qui 
Ol 
en 
01 
Ol 


O2 


02 
02 
01 
Ош 
Ou 
Qui 
O 
Ol 
On 


02 


№ = 12 


.68125674969 3980100 
.2024220944459 3580р 
с ОББ ооо ID 
-115315027696220 20D 
.629471055940404600 
57114093926 Олор 
.271149392612676500 
.62947105594040800р 
of (53250215962 4060D 
» 110776654 36970540D 
.20242269444595060D 
.6951256749594093100 


N= 13 


.79350759008819390р 
.23443251131591050р 
.12694451159465570р 
.574368882746129500 
-694179329375204 30D 
611039560984 34 350D 
58652 343749999S90D 
611034560984 34460D 
694179329 37626740D 
.87436888274613890р 
.12694451159465950р 
.23443251131591910р 
.7935075906884 3100р 


75 


-0. 


Ü 
-Q 


-0 


: D3RIVATIVES (continued) 


N = 14 


.91439461254679590D 
.265916927619570000 
.144 3976918035 35000 
.982087627944 362600 
.76617260993701340р 
.6586 3903475806550р 
613219721199 0 50D 
.61321072119930300D 
.658539034 75 8С66101 
„760172609937018900 
.922087027944 376900 
.144397691803540500 
.006591092761959250D 
.9143946125471 3950) 


N= 15 


10439175433549440D i 
. 30587 36445766 2590D 
a 
.10982495249623420D 
.84491321951164820D 
.71328504764050950р 
zt 
.62841796874999980D 
26. 
,71328584764051460П 
SON 
. 10982495 24962 36 /OD 
. 103127994164507 30D 
0. 30587 3644576659400 


16 312799416455820D 
64820514551937150D 
64820514551937180D 


Ола | ИІ) 


104 39175433554400D 


02 
02 


OL 
Ol 
Cd 
Ol 
Ol 
Ol 
OT 
Ol 
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APPENDIX F 


Coefficients of the Polynomials 
P i | 
y 09 / Geox) 
x 
The coefficients of the polynomials Ры (x) / (xx) are 
tabulated in order of smallest to largest root, X, and lowest to 
highest power. For example, 
х 1 
Р. (х) /(х-0.5000) = -1.05х°+30.4х'-188.5х2+389.8х2-252х*. 
These are the coefficients of the polynomials 


* 
P. (x) N-1 
N e K 
x-x, ` pr Я 
i K=O 


х 
where X. is a root of the shifted Legendre polynomial PN (x), and 


the coefficients were obtained by synthetic division 


5i.N-1 ^ CNN 
DEN SUE E-1US4 P. N-1 
NEN run 
b = C_.+x.b Й 
1350 МІ 111 , і- ІЛЕС 


The Legendre coefficients C,,, are defined by 
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* N M 
Р. (х) = > С : 
N m NM 
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PO (x)/(x-x,) (continued) 
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422000257 38625690D 
Ли 3[(55elp22490 NOD 


СОЕЕТТСТЕНТ5 ОЁ ТНЕ POLYNOMIALS 
* . 
P, (x)/(x-x,) (continued) 


(N = 12) 
I = 6' 
-0.17774090110084440D 


0.274116622955v0000D 
-0.10187901165549110D 
0.15981057468493100D 
-0.13172190226655770р 
0.63696597527710040D 
-0.1916664785 3966690D 
0. 3686461491 3043360D 
-0.45 307106764057 200D 
0.344061829119526102 
-0.1470353266 3486610D 
0.27041550000000000р 


AL 7 


-0.22563241685604960р 
O. 3514 3929209705040D 
-0.12928158809260100D 
0.199303087335 32540D 
-0.10040779750371220р 
0.75375947230744230р 
-0.21984244415561060D 
0.40977906012453700D 
-0.48875 309586617870р 
0.360994 36277086250р 
-0.15042183336513350D 
С.270415600000000001 


I= 8 


-0.316371346694904'70D 
0.49835 3021794298810D 
=0,17471512020309010D 
0.26141286017492700D 
-0.20231540762048540D 
0.9104 3508787125C80D 
-0.254640498807 905 80D 
Ome die 9225 3791910D 
-0.52710180022588890D 
O. 37820 359472091660D 
-0.15370194877502440D 
0.27041560000000000р 


OL 


07 


86 


Ер 
.73254291300179080р 
зад 
.35126118341389720р 
= 
.110769032932862500 
-0. 
.5069166 3294906840D 
o 
. 39445899926225400D 
.15666957684995930D 
.27041560000000000р 


-0 


-О. 
12803975 480242770D 
E 
513045 38720347610D 
-0. 
.13295684596250260р 
-0. 
„5545624 32956037500 
.59965118501394010D 
0 - 4084777154 3197 920D 
2704156 0000000000D 


-0 


=O. 
251858445685 36790D 
.660479439208806720р 
.70128491360441830р 
-0. 
„153 1633171 ГОД Да ВО 
aon 
- 59351356307 366690) 
-0. 
.41907005475455270р 
-0. 
.21041560000000000D 


-0 


(М = 12) 


I9 
4846 3460447902640D 


25556997895957640D 
2615276444 3786540D 
29447526 345421590D 
56525292661719900D 


I = 10 
865197411572201600 


41074814193766260D 
3312222557700330 
3352720859225 (50) 


159138264579192800 


I = 11) 
20858810712881270р 


416 37333025605900D 
37061910459839990D 
62655520668857870D 
1609529505175 8860D 


Ol 


COEBEPSICIENTS OF 
P(x) /(x-x 


(N = 12) 
I = 12 


-0.10846 360047086500D 
05155969046 35 14850D 
-0.92197417744400200D 
0.957142524459519300 
~0.47460933446 364170D 
O.1679045 3 36 3154500D 
-С.393378046172721800 
0.61/625 7799175 16800 
-0.64272493552206250р 
0.42529277095545140D 
-0.162000045 39231030D 
0.270415600000000000 


N = 13 
I = 1 


-0.10079715286689510D 
0.192434809108351200 
-0.907139761206638700 
0.153300979979990200 
-0.156074210727458600 
0.95417344 384 7699007 
-0.373118573646653800 
0.96475435951292820р 
-0.166734811715411700 
0.19042647081912740D 
-0.137866720 36695810D 
0.57285552860546860D 
-0.104005999999999801 
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D 


THE POLYNOMIALS 
.) (continued) 


(N = 13) 


I = 2 


—-—0.10429712616983220D 
0.18873297851273380D 
-0.83450914693393850D 
0.15833858496354280D 
-0.16096815846966000D 
0.98220090495688290D 
-0.38315116413652590D 
O.98778976478833180D 
-0.17012180818179740D 
0.19350609793310380D 
-0.13944189578075440D 
0.576 316130445 20726) 
-0.10400599999999980D 


I = 3 


-0.11101378477288740р 
0.20081267981965990р 
~0.88690991095794720D 
0.16795 375905 DD IIO 
-0.17026697 3696 30140D 
O.10351309566871580D 


-0.40194985295797750р | 


0.10305 335658962990D 
-0.176 33237685047210D 
0.19907 316154864720D 
-0.142241643712706 30D 
0.5823515 34547595 10D 
-0.10400599999999980р 


ngu 


-0.121776771004080401 
0.220150765672745 20р 
-0.970542507971309500 
OSEDS2TOT2359925059] 9E 
-0.18491574719641470D 
0.11176499795303420D 
-0.43087072104590820D 
0.1095 226767 350287 CD 
-0.185549991838031901 
0.20714875218926850D 
~0.14619803585958360D 
0.5906 31907 30108 370D 
-0.10400599999999980D 


COEFFICIENTS OF THE POLYNOMIALS 
* 
ps (х)/(х-х.) (continued) 


(N x 13) 


І - 5 


-0.13807456076844240D 
0.24938924169703470р 
-0.1096 396 34062665106 
0.20600177182381960р 
-0.20651787550969150р 
0.123740173894423201 
-0.47200090695790450р 
0.11850142742212430р 
-0.197985 3920225 3470р 
0.21770203249501490D 
-0.15119280842328430р 
0.60071303146785120р 
-0.1040599999999980р 


I = 6 


-0.162541060835 38880р 
0.29318277109067540D 
-0.128355'70496802190D 
0.239462090597264 30D 
-0.23767360347638980р 
0.14060256521493990D 
-0.529230087230211200 
0.130360062160615300 
-0.213792805711234401 
0.23058480315709960р 
-0.15704404056613840р 
0.612051476195 38480р 
-0.10400599999999980р 


I-7 


-0.20000000000000000D 
0. 35999999999999990D 
-0.15659999999999980D 
0.28899999999999990р 
-0.2825 3999999999980р 
0.16403231999999980р 
-0.60310559999999990D 
0.14542751999999990р 
-0.232928219999999801 
0.24545415999999970D 
-0.16349743199999980D 
0.62403599999999980D 
-0.10400599999999980D 


88 


(М - 13) 

I = 8 
-0.25989495 325629240D 
0.46625427624523020р 
-0.20017362 356530651006 
0. 36407741793960490р 
-0. 3476413554 32874800 
0.196236 34618423160р 
-0.70002386979225560р 
0.16379082 3362598000 
-0.25495904554972060р 
0.26170448342966820D 
ООС ПОСИ оа Шиит) 
0.6 3602052360461440D 
-0.10400599999999960D 


I = 9 


-0.362642558872232600 
0.64605549459456770D 
-0.27354641374456070р 
0.48160875083736360р 
-0.448457012652345700 
0.2394 327866 3642210р 
-0.920124609655950800 
0.18499071129263050D 
-0.27887545939564370D 
0.278436 37866464460D 
-0.17684807380864750р 
0.6473589685 3214790р 
-0.10400599999999980р 


І = 10 


-0.55920409476264 360D 
0.98648189703509230р 
-0.40282425526381880р 
0.67036126249861210D 
-0.576729607016894390D 
0.29412607766919960р 
-0.958807 35722672 390D 
0 .20770874136904660D 
-0.302993365 34 8000101 
0.29448811625060740D 
=0..10294253 (e203 9g oD 
0.657440092699915400 
-0.104005999999999200 


01 


COEFFICIENTS OF THE POLYNOMIALS 
P. (x)/(x-x;) (continued) 


(N 13) 


I = 11 


-0.10079532080864670D 
0.17328778717481680D 
-0.65084769641774720D 
0.95831383449131930D 
-0.749298470371563507 
Ода 1 8256 //5500 
-0.11024115446467420D 
0.22964 364791696 370D 
-0.32504456994910890D 
0. 3085418 3556664090D 
-0.18809455571035100D 
0.66572046545240380D 
-0.10400599999999980D 


I = 12 


-0.24271373142518990р 
0. 3828290 35771508000 
-0.105864682267220400 
0.13178219165 5391701 
-0.931724820382678700 
0.41498412167982510D 
-0.12281557116129330D 
0.247780 352 38250510D 
-0.34250008457189330D 
0.31929514242580010D 
-0.1919315522910261059 
0.67175 386955479200 
-0.10400599999999980D 


I = 13 


-0.12644666618088360D 
0.70245 338566609130D 
-0.147369 308371916 30D 
0.1617 3403145188670D 
-0.10667000989156780D 
0.455226978961970001 
—0.13097524020902070D 
0.259099 308307935100 
-0.35307976886615020D 
0. 32566732552037570D 
-0.19416523890094200D 
0.6752164713945 30 30D 
—0.10400599999999930D 


89 


-0. 
,2104 362867045 3990D 
-0. 
.238371254753723300 
.20441945032009650D 
„2057560 14 ЗОРО 
-0. 
. 30468070295618890D 
-0. 
.982637659099974607 
. 987 35094411646760D 
6 395 2860422 322550D 
2409747 23480016 60D 
.40116599999999990D 


-0 


-0 
-0 


-0. 
.216717599647022 30D 
-0. 
.24519315203176540D 
-0. 
.211122677502515901 
98820046 304836740D 
. 3113954 3910837490D 
-0. 
. 9985600121101071 30D 
-0. 
.64537983781745240р 
-0. 
.40116599999999990D 


-0 


N= 14 


I = 1 
1006 9054782458 360D 


10783518111968030D 


96481786447879590D 
66011574 395699110D 


I = 2 
10371104684228440D 


11100486193258090D 
29224391985046560D 


67295 358233064620D 
100020710 37698320D 
24213507457 301090D 


COEFFICIENTS OF THE POLYNOMIALS 
Р. (х)/(х-х.) (continued) 


'(N = 14) 


EIL 


~0.10945701042210040D 
0.22866164573420590р 
-0.117024197236 39010D 
O" 25 JI OO J 5571141]0D 
-0.30702047934305 320D 
0.22120224016802250р 
-0.10319175563517800р 
О 82 Бета пе злбелот 
-0.6965916 3832994400D 
0.10282640822576100р 
-0.10233718266634350D 
0.65578080356590450р 
-0.24410564786193610р 
Ооа 25000 2999990р 


I = 1 


-0.11853306384700660D 
0.24751442830826710D 
-0.126504242876347007 
On 270390 1Or 4 ome 3 20D 
-0.33007450983788320D 
0.23661141365139410D 
-0.10989809578462980D 
0.34276783051801440р 
-0.73194343166031770р 
Q, TO (162360 719067210D 
~0 .10568339927924010D 
0.67049367136484310D 
-0.24697197272741550D 
0.40116599999999990D 


90 


-0. 
.21544005597036710D 
-0. 
. 30016075444191710D 
. 30862 339215554540D 
.259271272140 060800 
-0. 
: 359129291956432900 
-0. 
SIN 20580) 229) ү Лр К) 
-0. 
.68912423001244720р 
-0. 
4011652999998009990D 


-0 


-0. 
.316097145128032102 
-16077330637 790D 
35090667 15111745702 
-0. 
.29060704866 359880D 
-0. 
.40429806385 8950401 
-0. 
ТООС Оа 3T 998 
-0. 
„711087935685351401 
-0. 
Os 


229, 


(N = 14) 


I = 5 
131991539024379801 


14049910526 392 390D 


11941162739159840D 
78025480169954270D 
1100396 3250961380D 
26042288827 7337200 


І = 6. 
1516171060 74675200 


411208838051084 30D 
1324156786 1503020р 
84294751415744180D 
11534911543210250D 


25435704837033370D 
40116599999999990D 


01 
03 
05 
06 
ОТ 
08 
Og 
09 
6/5 
10 
1) 
09 


08 


COEFFICIENTS OR THE PODYNOSATS 
Р “(х)/(х-х.) (continued) 


(N = 1!) 


I = 7 


-0.180496464035 37720р 
оо 1846 1614025 420) 
-0.190319357998493001 
Q.4124] 29 OP ии веш 
-0.47843243277376570D 
0. 33382697661047 380D 
-0.1498 36580642135 30D 
0.44986 2851354082 80D 
-0.92126560909521140D 
0.12900962569365 4401 
-0.12149596007792130D 
0.73555445760512930D 
-0.25859050142234330D 
Ompuguso599909999990D 


I = 8 


-0.22422905954241500D 
0.416505 315199012260) 
-0.23441235145503930D 
0.502449707039155000 
-0.5'1416 32379045 8560D 
0. 3932444851067 3690D 
-0.17281744290778780р 
0.50733061810099180р 
-0.101562016907 338300 
O.139110410451295 20D 
=0.1282843939002538020 
©. 7615932405 3/00740D 
-0.26292529857765640р 
САО аззз999в0й 
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-0 


-0 


-0 


(не 14) 
I = 9 


.2937 342240839 3690D 
.608213891132445501 
-0. 
6 3808107 291 382380D 
-0. 
4145647 3659044190D 
-0. 
„577402630890027501 
. 112461046915 IT OD 
„15О216030945 06.800 
-0. 
.727899155241347100 
-0. 
011659999 II Б 


3028924491587 7940D 
71229930356286800D 
2024998052 3924240D 


1354258545 T оле 
20715875162966600D 


I = 10 


.412582645456995807 
"8494011115 2544690D 
70164919835908124590D 
7058311048055 315070D 
7016581145 1818192800D 
т 
70:813144370684400100 
70:40116599999999990р 


41549543307351450D 
91315 3342897659700 
23940203695940500р 
ра доиде а OD 
1425407179 9220D 
27109291172266260D 


COEPRICIENBSCOPF THE POLYNOMIALS 
+ 2% ; 
Бы (x)/(x-xi) (continued) 


(N= 14) 
I ж. 11 


-0.639576149941905000 
0.130220414972016 30D 
-0.6151312840847 3890D 
O.11896555879521140D 
-0.11981847811414090D 
QE EOD QI 
-0.28231372004949410D 
0.746764251576655500 
-0.13655274325948350р 
0.1729970004 30881200 
-0.1491793752215510D 
0.835924 7986 3585 380D 
-0.27454382727258420р 
0.40116599999999990D 


I = 12 


-0.115741621575861500 
0.22966128234427870р 
-0.99808481529113370D 
0.17096 37067 3251690D 
-0.156647 364905725 30D 
O.88131683554205400D 
-O. 3271390225823 3050D 
0.832575056 30568550D 
-0.14779539988498060р 
0.18294164821220520р 
2015406310536 T3712 TOD 
0.854887829222670900 
-0.27735015213806370р 
0.40116599999999990D 


Ol 
04 


09 


32 


(N = 14) 


I = 13 


-0.27946576521975180р 
0.50877699303161640р 
-0.16299086 393618340] 
0.236229996934389300 
-0.195838438202724901 
0.10327411926099290D 
-0.366719469226 31620D 
0.904169 30565976880D 
-0.15678030641774040р 
0.19062727573007680р 
-O. 19904031902 0955 2p 
0.86885 37429413 3100D 
-0.27938072542698890р 
0.40116599999999990D 


І = 19 


-0.14581307272487430D 
0.93592930947540420D 
~0.22757146947683940D 
0.29087565 319279450D 
-0.22499462811690370D 
0.11370404087734690D 
—-0.39256 36589841 3480D 
0.94914118241443470D 
-0.16226264769256810D 
0.19521312846090240D 
-0.16164795526579080р 
0.87592612246302460D 
-0.28054107651998310D 
0.40116599999999990D 


COEFFICIENTS OF THE POLYNOMIALS 
P.C (x)/(x-x,) (continued) 


-0.100604015247729700 
0.240437 40422705070 
-0.14124 3615240087 80D 
0.359312859820643200 
-0.49612129778936520% 
0.41848546058240950D 
-0.23113042032021540D 
O.Ə7155279695345 eo) 
-0.22974108527756610р 
0.42722964284119880D 


-0.5577274 3622680000D . 


0.499894641722260740D 
-0 . 292829058096 82000D 
0100919516541309 30D 
OTTS 11753T99999J990D 


I=2 


-0 . 1032 37 81804559840D 
0.246705109105339801 
-0.144876762835258107 
0.36834403466923620D 
-0.50817764644264680р 
0.42619702805791220D 
-0.236178969898214207 
0.88915862209255000D 
-0.23393827394690690р 
0.43408171474656100D 
-0.56525518329132970D 
0.50521073445626440D 
-0.29500729908815030р 
0.10131288779044060D 
EUSb IP953:999999990D 
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(N 2 15) 
I = 3. 


-0.10821295604528690D 
0.258540222271777400 
—0.1517 304 38616 39020D 
. 385354060887 77140D 
-0.530826584 33782880D 
.44638133558104910D 
-0.24559287141861270D 
921820937 36457290D 
-0.24161749132397410D 
446625679597 389 30D 
-0.57892115871097 7001 
.51476895013676600D 
-0.29888074454704940D 
0.10200367891668490D 
-0.15511751999999990D 


I =L 


-0.11598116781339040D 
0.27700966416469860D 
-0.16240833159132670D 
0.41177960837793960D 
~0.5658676 90049055 30D 
0.47436110584700210р 
-0.25997785262680260р 
0.97131080496675650D 
-0.25328663064844250D 
0.46522268248212260D 
-0.59891050309551950D 
0.5285 3821741328200D 
-0.30436640474547060D 
0.10296 376990336900D 
-0.15511751999999990D 


ev (ev muere eed 


01 
03 


06 
OT 
08 
09 
09 
10 
10 
10 
1110) 
110 
10 
09 


COBEBSIEC TI ENDSESORP THE POLYNOMIALS 
* ` | 
В. (х)/(х-х.) (continued) 


(N = 15) 
T= 5 


-0.127 3097077 3184230D 
0.30392251000383360р 
-0.177929026795602200 
0.450023 37213037700D 
-0.61627019399481760D 
0.51427785627769750р 
-0.28028954072007160р 
0.1040325 3648001450р 
-0.269240316547203200 
0.49035297017800650D 
-0.62541611670470500D 
0.5464192565 3287 200D 
-0.31133126064361470D 
0.10415 357462979740D 
zO opi (00999990D 


I = 6 


~0 1434564527 3585880D 
0.342237499497777001 
-0.19994618991256290D 
0.50394152936850410D 
-0.68669710511210130р 
0.56940074737338770D 
20.30793105169337000D 
0.11326167563691360D 
-0.290148 38985 33506 30D 
0.522550674616574400 
-0.65054023344066 340D 
0.56619159222852080D 
-0. 319583577835 30470D 
ОТ 225 275026933400 
-О du/51999999990D 
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(N = 15) 


I= 7 


-0.1605009860 3607 240D 
O. 3968 3010652589540D 
-0.231156146418811000 
Oro TOCIT a Ol 1 
-0.78439327865484000D 
0.64460155784839010р 
-0.34487542895894950р 
0.125 30451403 3459201 
-0.31670575123719550D 
0.961227 301545203220) 
-0.69819660060510230D 
0.59344986 9885 32170D 
-0.32886568640604680р 
0.10702182755502340р 
-0.15511751999999990D 


I = 8 


-0.20000000000000000D 
0.476000000000000001 
-0.27607999999999990D 
0.687 34400000000000D 
-0.92067919999999990D 
0.74703439999999970D 
-0.39377575999999980D 
0.1407 3497599999980D 
-0.34956416599999970D 
0.60968317199999950D 
-0.743850911999999601 
0.6215 3985599999960D 
-0.33685154799999980D 
0.10858226399999980D 
-0.15511751999999990D 


01 
O3 
05 
06 
07 
09 
09 
10 
10 
10 
10 
10 
10 
1O 
09 
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ELEMENTS OF THE EXPLICIT INVERSE 
MATRIX WITH DIVISION BY WEIGHTS 
From equations (37) and (45) it is evident that the explicit 
inverse matrix with division by weights is the necessary matrix for 


use in obtaining the temperature function: 
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